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ABSTRACT 


The initial success o f the Rigidly Rotating Magnetospher e (RRM) model application to 
the B2Vp star a Ori E by [Townsend, Owocki & Gro^ ( 2005| ) triggered a renewed era of ob¬ 
servational monitoring of this archetypal object. We utilize high-resolution spectropolarimetry 
and the magnetic Doppler imaging (MDI) technique to simultaneously determine the mag¬ 
netic configuration, which is predominately dipolar, with a polar strength Bd = 7.3 — 7.8 kG 
and a smaller non-axisymmetric quadrupolar contribution, as well as the surface distribution 
of abundance of He, Fe, C, and Si. We describe a revised RRM model that now accepts an 
arbitrary surface magnetic field configuration, with the field topology from the MDI models 
used as input. The resulting synthetic Ha emission and broadband photometric observations 
generally agree with observations, however, several features are poorly fit. To explore the pos¬ 
sibility of a photospheric contribution to the observed photometric variability, the MDI abun¬ 
dance maps were used to compute a synthetic photospheric light curve to determine the effect 
of the surface inhomogeneities. Including the computed photospheric brightness modulation 
fails to improve the agreement between the observed and computed photometry. We conclude 
that the discrepancies cannot be explained as an effect of inhomogeneous surface abundance. 
Analysis of the UV light variability shows good agreement between observed variability and 
computed light curves, supporting the accuracy of the photospheric light variation calcula¬ 
tion. We thus conclude that significant additional physics is necessary for the RRM model 
to acceptably reproduce observations of not only a Ori E, but also other similar stars with 
significant stellar wind-magnetic field interactions. 


Key words: stars: magnetic fields - stars: rotation - stars: early-type - stars: circumstellar 
matter - stars: individual: HD 37479 - techniques: spectroscopic 
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National Research Council of Canada, the Institut National des Sciences de 
TUnivers of the Centre National de la Recherche Scientihque of France, 
and the University of Hawaii 
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1 INTRODUCTION 

A small percentage of main sequence A- and B-type (Ap/Bp) stars 
are characterized by peculiar surface abundances of chemical ele¬ 
ments and strong (0.1-30 kG), large-scale magnetic fields. These 
magnetic Ap/Bp stars exhibit periodic variability of spectral lines, 
photometric brightness, and longitudinal magnetic field strength 
that is explained in the context of the Oblique Rotator Model 
(ORM) described by |Stibbs| ( |1950| ). The peculiar element abun- 
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Figure 1. Comparison between observed (black) and computed (red) Stokes 
/ and V profiles of the He i 5876A and 6678A lines. Rotational phases are 
indicated to the right of the Stokes V profiles. Spectra corresponding to 
consecutive phases are shifted vertically. The bars at the lower left corner 
of each panel indicate the horizontal (1 A) and vertical (5 per cent) scales 
of the profile plots. 


dances present themselves in the form of inhomogeneous, asym¬ 
metric “spots” on the surface of the star, causing spectral line shape 
and strength variability, and in some cases photometric brightness 
variations, all modulated on the stellar rotation period. These inho¬ 
mogeneities on the surface of magnetic chemically peculiar stars 
are thought to be formed by direct or indirect interactions between 
chemical diffusion and the stabilizing influence of magnetic field 
lines (see e.g., |Michaud|1970| ). 

The technique of Doppler imaging (DI) was originally used 
by |Goncharskij et al.| ( |1982| ) to map the abundance of a chemi¬ 
cal element on the surface of a star assuming rotational variabil¬ 
ity of line profiles with time. This simple, yet elegant method of 
studying stellar spots was eventually extended to simultaneously 
map both chemical abundances as well as a star’s surface magnetic 
field. This method, magnetic Doppler imaging (MDI), was previ¬ 
ously employed to study Ap/Bp stars, and successfully investigates 


both the surface inhomogeneities of various elements and the mag¬ 
netic field topology. The MDI code iNVERSlO was developed by 
[Piskunov & Kochukhov| ( |2002| ) to accurately model stellar mag¬ 
netic fields and surface abundance features from an inversion of 
high resolution spectropolarimetry. The code was originally devel¬ 
oped to operate with data in all four Stokes parameters. However, it 
is also able to handle a set of data comprised of only Stokes / and 
V spectra. 

Using abundance maps derived from DI, |Krticka et al.| { [2007| ) 
showed that both the photometric and spectroscopic variability of 
another He-strong star, HD 37776, can be explained by the pres¬ 
ence of inhomogeneous abundances across the stellar surface. Sim¬ 
ilar work corroborates this claim for the case of the rapidly rotating 
Ap star CU Vir ( jKrticka et al.||2012| ). [Krticka et al.jpOT^ deter¬ 
mined that the photospheric spots of HD 64740 are not sufficient 
to cause detectable photometric variability. Further, these authors 
identified UV line variability due to these spots, without any evi¬ 
dence of circumstellar features. 

a Orionis E (HD 37479), the prototypical magnetic He-strong 
Bp star, has long been known to exhibit variable line profiles, indi¬ 
cating an inhomogeneous surface distribution of various elements. 
The discovery of its magnetic field by [Landstreet & Borr^ ( |1978j ) 
suggested that this star was a higher mass extension of the Ap/Bp 
phenomena, and could be physically represented by the ORM. A 
study by [Reiners et al.jpOOOj ) modeled historical and new spectra 
to investigate line profile variations to determine their relationship 
to local differences of element abundances. The authors found re¬ 
gions of He overabundance and metal (C ll and Si III) underabun¬ 
dance slightly offset from each other, suggested to be located at the 
magnetic poles. 

[Walborn[ ( [1974[ ) had discovered broad, variable Ka emission 
features, which, when considered in the context of the ORM, in¬ 
dicated that the emitting material was magneto spheric in nature 
( [Landstreet & Borra|1978[ ). Various authors worked to develop re¬ 
alistic magnetospheric models for this and other similar stars, with 
the most recent being the development of the Rigidly Rotating 
Magnetosphere (RRM) model by [Townsend & Owocki[ ( [2005] l. The 
model analytically describes a rapidly rotating star surrounded by 
plasma trapped in its magnetosphere. Although the RRM model 
was relatively successful in reproducing the magnetospheric spec¬ 
tral and photometric variability of a Ori E ( [Townsend, Owocki &[ 
[Groote|2005[ ), major inconsistencies remained. 

To begin addressing these issues and the revision of the RRM 
model, [Oksala et al.[ ( [2012[ hereafter referred to as Paper I) pre¬ 
sented and analyzed high resolution intensity and circular polar¬ 
ization spectra a Ori E. These spectra identified various spectral 
lines that exhibited significant variability and presented an updated 
magnetic field characterization, suggesting a more complex, higher 
order field structure. This result was in strong contradiction with 
the offset dipole configuration adopted for the specific case of a 
Ori E by [Townsend, Owocki & Groote| ( [2005[ ). 

A remaining issue with the RRM model, which we aim to ad¬ 
dress here, are persistent discrepancies between the observed opti¬ 
cal photometric light curve and the model light curve (see for exam¬ 
ple Eigure 8, this paper). As noted by [Townsend, Owocki & GrooteJ 
( [2005[ ), a significant observed feature of increased brightness is un¬ 
matched at rotational phase 0.6 (see their Fig. 1). As these authors 
invoked the de-centered dipole as a means to mimic the unequal 
depths of the light minima, a change of the magnetic configura¬ 
tion would cause changes to the synthetic light curve, and further 
differences will likely arise in the model-observation comparison. 
The most precise photometric data to date were obtained over a pe- 
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Figure 2. Surface magnetic field distribution derived for a Ori E with magnetic DL The star is displayed at five equidistant rotational phases as indicated at 
the top of the figure. The rotational axis is indicated by the short vertical bar, while the thick line shows the stellar rotational equator. The maps correspond 
to an inclination angle i = 75°and a Urot sin i = 140 km s“^. a): Maps of the magnetic field modulus. The contour lines are plotted with a step of 1.0 kG. 
b): Maps of the radial magnetic field strength. The contour lines are also plotted with a step of 1.0 kG. c); Orientation of magnetic field vectors. In the vector 
maps, light (red) arrows point outward from the stellar surface. Dark (blue) arrows correspond to vectors pointing inwards. The arrow length is proportional to 
the field strength. 


riod of three weeks with the MOST micro-satellite, and were pre¬ 
sented and analyzed by [Townsend et al.| l |2013| ). The emission and 
absorption features previously seen in the light curve of a Ori E are 
clearly present, and support the notion that these features are real, 
and also stable. The properties of the two minima strongly differ, 
with one eclipse deeper than the other, while the more shallow fea¬ 
ture is present during a longer portion of the rotation period. The 
excess brightness at phase 0.6 is quite evident, and the high preci¬ 
sion reveals a gradual decline of this excess before returning to the 
“continuum” level. 

This subsequent work continues our re-evaluation of the RRM 
model, utilizing cutting edge techniques to analyze the magnetic 
and abundance structure of a Ori E via MDI, revise the RRM 
model, and, through a thorough light curve analysis, address any 
remaining discrepancies between model and observation. Section 2 
explains the MDI procedure and results. Section 3 presents the re¬ 
vised RRM model. We discuss in Section 4 the light curve analysis 
of the optical photometry. An analysis of UV variability is explored 
in Section 5. Section 6 presents a summary of this paper. 


2 MAGNETIC FIELD AND ELEMENTAL ABUNDANCE 
MAPPING 

2.1 Magnetic Doppler imaging (MDI) 

Eor input into MDI, we obtained a total of 18 high-resolution 
(R=65000) broadband (370-1040 nm) Stokes / and V spectra of 
a Ori E. Sixteen spectra were obtained in November 2007 with 
the Narval spectropolarimeter on the 2.2m Bernard Lyot telescope 
(TBL) at the Pic du Midi Observatory in Erance. The remaining 
two spectra were obtained in Eebruary 2009 with the spectropo¬ 
larimeter ESPaDOnS on the 3.6-m Canada-Erance-Hawaii Tele¬ 



Figure 3. Longitudinal magnetic field measurements (red filled squares) 
for (7 Ori E (Paper I), as well as data reported in Landstreet & Borra (1978) 
and Bohlender et al. (1987) (open circles). Each data point is plotted with 
corresponding Icr error bars. The solid curve is the longitudinal magnetic 
field curve computed using the MDI-derived magnetic field structure. The 
data are phased using the ephemeris of Townsend et al. (2010). 


scope (CEHT), as part of the Magnetism in Massive Stars (MiMeS) 
Large Program ( jWade et al.|2011| ). The observations and their re¬ 
duction are described in more detail in Paper 1. 

In this work, we utilize the iNVERSlO MDI code 
Kochukhov|2002 1 and its recent derivative invers13 
Wade & Shulyak 2012| [Kochukhov et al.|2013| l to reconstruct the 
surface abundance distribution and magnetic field topology of a 
Ori E using these high-resolution Stokes / and V spectra. Both 


Piskunov & 
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Figure 4. Comparison between observed (black) and computed (red) Stokes / profiles of lines used for chemical abundance mapping.. Rotational phases are 
indicated to the right of the profiles. Spectra corresponding to consecutive phases are shifted vertically. The bars at the lower left corner of each panel indicate 
the horizontal (1 A) and vertical (5 per cent) scales of the profile plots. 


codes use least-squares minimization together with a regulariza¬ 
tion procedure to fit synthetic spectra to spectropolarimetric obser¬ 
vations. The model spectra are adjusted until the solution converges 
to the simplest surface distribution that properly describes the ob¬ 
servational data. 

Because only Stokes / and V spectra were obtained, the final 
magnetic field solution is not unique, and the choice of regulariza¬ 
tion becomes important. As described by [Piskunov & Kochukhov| 
( |2002| ), without Q and U spectra, we must use an additional a priori 


constraint, multipolar regularization, which constrains the results 
according to agreement with a best fit multipolar model magnetic 
field. Additionally, Stokes IV observations only provide informa¬ 
tion about the global magnetic field geometry, and is thus unable to 
reveal small-scale magnetic field structures, such as those seen in 
some full Stokes IQUV inversions of Ap stars (e.g.,|Kochukhov et 
al.|20041|Kochukhov & Wade|20 101 [Silvester, Kochukhov & Wade 
20\4\ We also include the longitudinal field curve as an additional 
constraint in magnetic inversion following the procedure described 
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Figure 5. The chemical abundance distributions of C, Fe, and Si derived from Stokes / line profiles. The star is shown at 5 equidistant rotational phases viewed 
at the inclination angle, i = 75° and v sin i = 140 km s“^. The scale gives abundance as CEiem corresponding to log(NEiem/Ntot) relative to solar values. 
The rotation axis is vertical. The contour step size is 0.5 dex. 
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Figure 6. The chemical abundance distribution of He derived from Stokes / line profiles of the 4713 A He i line. The star is shown at 5 equidistant rotational 
phases viewed at the inclination angle, i = 75° and vsini = 140 km s“^. The scale gives abundance as log(He/H) relative to solar values. The rotation 
axis is vertical. The contour step size is 0.5 dex. 


by |Kochukhov et al.| \2002\ . The resulting field configuration is 
computed as part of the multipolar regularization procedure. Thus, 
the multipolar parameters reported here represent an approximation 
of the 2-D MDI magnetic map, which is derived without using the 
spherical harmonic parameterization. 

2.2 Magnetic field topology 

To determine the magnetic field topology of a Ori E, the line pro¬ 
files of the He I lines at 5876 A and 6678 A were chosen for the 
inversion process. These lines both have strong magnetic signa¬ 
tures in their Stokes V spectrum, while other possible lines exhibit 
much weaker polarization features, and are therefore less suitable 
for use in magnetic inversion. We use as input the parameters (Tgff, 
log g, and Urot sin i) determined in Paper I, which are listed in Ta¬ 
bled The rotational period, Prot, was derived by [Townsend et al.| 
( |2010| ). While these parameters are rather well established for a Ori 
E, there is a large uncertainty in the determination of the inclination 
angle, i. This is a result of uncertainty in the size and distance of 
the star, and as such the radius. A range of radius values have been 
reported via various methods, and so we compute MDI maps at 4 
fixed inclination angles to understand the effect of this parameter. 


Corresponding to a range in radii of 3.3-4.0 R*, models were com¬ 
puted for i = 55, 65, 75, and 85°. This range is adopted based on 
the discussion of the stellar radius in the Appendix of |Townsend et| 
[^ ( |2013| ), accounting for approximate uncertainties in the angular 
diameter (|Groote & Hunger] 1982|l and the cluster distance (|Sherry| 
let al.|200^ 


While the inversIO code relies on LTE model atmospheres, 
this assumption is not always realistic. Eor hot stars, the LTE as¬ 
sumption begins to break down starting at Tgff > 15000 K ( |Mi-| 
|halas & Athay|197^ , beginning with visibly different line profile 
strengths and/or shapes. These effects become particularly notice¬ 
able above Tgff ~ 22000 K, resulting in growing discrepancies 
with LTE model spectra in Hi lines and He I lines, particularly 
those in the red-ward part of the optical spectrum (i.e., 5876 and 
6678 A; Przybilla, Nieva & Butler|2oTT] ). Less obvious effects are 
seen in metal lines and blue-ward He I lines; the study by |Przy-| 
|billa, Nieva & Butler| ( |201 1| ) finds little difference between LTE and 
NLTE models for these lines. Thus, to properly model the He I lines 
we have chosen, these changes need to be incorporated into the syn¬ 
thetic spectrum. The implementation of full NLTE radiative trans¬ 
fer into the INVERS suite of codes is not a trivial task, and we pro¬ 
ceed to recover the line profile shapes only for lines to be included 
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Figure 7. Visualizations of the plasma distribution predicted by the ARRM model, at the same 5 rotational phases adopted in previous figures. The star is 
colored according to the radial magnetic field strength, as in the middle panels of Fig.[^ selected field lines are plotted in green, and the magenta arrow shows 
the star’s rotation axis. The surrounding distribution of magnetospheric material is shown via an orange-colored volume rendering, with the more opaque parts 
corresponding to higher column densities and vice-versa. The tick marks on the surrounding boxes are spaced by 1 i?*. 


in the magnetic field determination. The effect of NLTE departure 
is more severe for red He I lines, however, the magnetic signature is 
also stronger in these features, making them useful for determining 
the magnetic configuration. As a solution, we have used the NLTE 
model atmosphere code TLUSTY ( |Lanz & Hubeny||2007| ) to de¬ 
termine the departure coefficients for the upper and lower atomic 
levels of the He transitions in question. The code was set to only 
consider particle conservation and statistical equilibrium to calcu¬ 
late both the population levels and their associated departure coeffi¬ 
cients. Departure coefficients were computed for various values of 
helium abundance with a LTE ATLAS9 ( [Castelli & Kurucz|2006| ) 
model atmosphere with Tgff = 23000 K and log^ = 4.0 as the 
initial input for the TLUSTY computation. Then, during inversion, 
iNVERSlO interpolated departure coefficients for the local value of 
He abundance and modified the line absorption coefficient and the 
source function accordingly. 

With the longitudinal magnetic field curve derived in Paper I 
as an additional constraint in addition to the multipolar regulariza¬ 
tion ( [Piskunov & Kochukhov|200^ , we obtained a good fit to both 
the Stokes I and V line profiles of the two He I lines, shown in 
Pigure[^ These figures correspond to an inclination of i = 75°. The 
best fit magnetic field configuration over the considered range of in¬ 
clination angles (z = 55 — 85°) suggests a magnetic field that can 
be approximated as a dipole with polar strength Sd = 7.3—7.8 kG, 


with obliquity /3d = 47 — 59° and phase angle 7 d = 97 — 98°, and 
a smaller non-axisymmetric quadrupole component with strength 
Bq = 3 — 5 kG. Figure shows the surface magnetic field dis¬ 
tribution maps derived for a Ori E corresponding to this topology 
(again, for i = 75°). The dipole is centered, but the quadrupole axis 
is misaligned with the dipolar axis, leading to a geometry where the 
positive and negative poles are clearly not separated by 180°. The 
asymmetry is further evident in the longitudinal field curve shown 
in Figure where the observed longitudinal field measurements 
from all epochs are plotted along with the synthetic curve computed 
from the MDI magnetic maps. The data are very well matched, and 
corroborate the conclusions of Paper I, particularly their Figure 7, 
which shows a comparatively poor fit to the observed data using the 
original RRM de-centered dipole field configuration. 


2.3 Metal abundance 

With a magnetic field topology corresponding to an inclination 
i = 75° set as a fixed parameter, inversIO was then used to com¬ 
pute the distribution of various chemical elements on the surface of 
the star. This inclination was chosen as it was used for the origi¬ 
nal (T Ori E RRM analysis by jTownsend, Owocki & Groote| ( |2005] ), 
and we expect far less uncertainty in the determination of the sur¬ 
face fields with varying inclination angle. The set of high resolution 
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Figure 8. Comparison between the observed (black dots; [Hesser et al.|1977t 
and modeled photometric light curves for i = 75° (red curve) and i = 85° 
(green curve), for the Stromgren n-band. 

time-series spectra of a Ori E contain lines suitable for modeling 
the abundances of Fe, Si, C, and He. Other chemical elements are 
present in the spectrum, but their line strength was too weak for 
accurate determination of any surface structure. He lines required a 
slightly more rigorous treatment, and will be discussed in the next 
section. 

The Stokes I spectral lines were fit, with the final results 
shown in Figure There is good agreement between the obser¬ 
vations and the model fits for all of these lines. Figure displays 
spherical maps of the abundance distributions of Fe, Si, and C. The 
maps of Si abundance were determined using three closely located 
Sim lines at 4552.62 A, 4567.84 A, and 4574.76 A. The surface 
distribution of C was computed from the C ll 4267 A triplet. Fe 
maps were produced using the Felll doublet lines at 5127.39 and 
5127.63 A and the Fe III singlet line at 5156.11 A. The surface maps 
shown in Figure [^indicate that Fe, Si, and C all show a similar pat¬ 
tern over the stellar surface. The minimum abundances are found at 
rotational phase 0.8, in both a spot at the equator and at the visible 
pole. An equatorial spot at ~ phase 0.6 has the maximum value on 
the surface. 


2.4 He abundance 

As discussed in Section \l?2\ the two helium lines used for mag¬ 
netic mapping are more sensitive to departure from NLTE. Further, 
these lines are very strong, and as such they have a weaker ability 
to diagnose changes in abundance. We find that He I 5876 A and 
6678 A also contain contamination from the circumstellar material, 
which was confirmed in Paper I by evaluating the variability of 
the equivalent widths (EWs) of these lines. We therefore choose 
to determine the photo spheric helium abundance using the weaker 
He I 4713 Aline, which does not show any distortion due to either 
NLTE effects or the magnetosphere. Further, instead of using the 


standard inversIO MDI code to compute synthetic He line pro¬ 
files, we chose to use a different imaging code, invers13 (for more 
details see [Kochukhov, Wade & Shulyak|2012l|Kochukhov et al.| 
|2013| ). This MDI code is capable of computing local Stokes param¬ 
eter profiles according to local atmospheric conditions, modified by 
the presence of chemical or temperature spots. In the case of a Ori 
E, which shows a large He overabundance and a non-uniform sur¬ 
face distribution of this element, we computed a grid of LTE model 
atmospheres with LLmodels code ( |Shulyak et al.|2004] | for a range 
of He/H ratios. This grid was then used in modeling of the He I 
4713 A line, allowing us to derive a more accurate He surface map 
than when using a single mean atmosphere. As with the metals, the 
He abundance distribution was computed with the fixed magnetic 
field geometry determined from the red He lines. 

Figure shows a He abundance map indicating a large spot 
of overabundance at phase 0.8, located in the lower hemisphere. 
This corresponds with the EW variation in the He I 6678 and 5876 
A lines derived in Paper I. The minimum abundance is located at 
phase 0.6, with a normal, solar level of He. The spot of enhanced 
He does not appear to be correlated with the location of the mag¬ 
netic poles. Also of note is the anti-correlation between the strength 
of metals and He, previously established in Paper I from EW vari¬ 
ations. 


3 THE REVISED RIGIDLY ROTATING 

MAGNETOSPHERE MODEL FOR cr ORI E 

The initial success of applying the Rigidly Rotating Magnetosphere 
(RRM) model to a Ori E, by [Townsend, Owocki & Groote| ( |2005| ), 
triggered a renewed era of observational monitoring of this archety¬ 
pal object, across a range of different wavelengths. In Paper I, 
we presented the first-ever high-resolution spectropolarimetry of 
a Ori E, which reveals that the star’s magnetic topology is more 
complex than assumed in the original RRM analysis. Here, we ex¬ 
plore the consequences of this revised topology, considering the 
analysis of Section |2^ for the RRM-predicted distribution of cir¬ 
cumstellar material. 

3.1 The Offset-Dipole RRM Model 

The RRM analysis by [Townsend, Owocki & Groote] ( [2005[ ) as¬ 
sumed a dipole magnetic topology tilted at an angle /3 = 55° with 
respect to the rotation axis and then offset by 0.3 i?* in a direction 
essentially perpendicular to both rotational and magnetic axes. The 
offset results in an mass/density asymmetry in the predicted mag- 
netospheric mass distribution (with respect to the plane containing 
both axes), which nicely explains the observed inequalities in the 
depths of the light-curve minima and the heights of the Ka emis¬ 
sion peaks. 

[Townsend, Owocki & GrooteJ ( [2005[ l reported that their offset- 
dipole topology produces a reasonable fit (reduced — 1-79) to 
the longitudinal magnetic field (Bi) measurements available at that 
time, only marginally worse than a centered dipole. However, a re¬ 
view of the field synthesis code used by these authors has revealed 
a couple of bugs, which tend to suppress the effects of the offset on 
the Bi curve. With these bugs fixed, the predicted Bi is clearly in¬ 
consistent with the observations: the transition from magnetic max¬ 
imum to minimum is gradual, while the transition from minimum 
of maximum is quite steep - precisely the opposite of the observed 
behavior (see Fig. 3). 

This mismatch could be partly mitigated by offsetting the 
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Figure 9. Comparison between observed (left) and modeled (middle and right) Ha dynamic spectra of a Ori E, in each case plotting the residual flux as a 
function of velocity and rotation phase. The middle plot is determined for an inclination angle, i = 75°, while the right plot is the same, but for i = 85°. The 
observational data and the synthetic reference spectrum used to calculate the residuals are described in Section 4 of paper I. 


dipole in the opposite direction (at the expense of worsening the 
light curve and Ho fits). However, recognizing that an offset dipole 
is at best a crude approximation to the star’s true magnetic topol¬ 
ogy, and that we now have access to a more physically realistic 
topology (Sec |2.2| ), a better approach is to discard the offset dipole 
approach and instead extend the RRM model to handle arbitrary 
non-dipolar fields. 


3.2 A New Arbitrary-Field RRM Model 


The Arbitrary-Field RRM (ARRM) model is an extension of the 
original RRM formalism to allow consideration of arbitrary cir- 
cumstellar field topologies. The theoretical basis remains the same: 
the distribution of density p along a given closed magnetic flux 
tube is calculated by solving the equation of hydrostatic equilib¬ 
rium subject to gravitational and centrifugal forces. This yields 


p(s) = 


po exp 
0 


r ^(s)-^o i 

[ ^ fcT J 


Sa < S < Sb 

Otherwise 


( 1 ) 


where s is the coordinate measuring arc distance along the flux 
tube; Sa and Sb are the coordinates of the maxima in the effec¬ 
tive potential T>(s) that are situated closest to the tube footpoints, 
and the other symbols have the same meaning as described by 
[Townsend & Owocki| p005jl. In this expressi on, which is equiv¬ 
alent to equation (25) of Townsend & Owocki ( [2005 j ), po = p{so) 
is the density at an arbitrary reference coordinate s = so, and 
T>o = ^(so) is the corresponding effective potential. 

To determine po, we integrate the density distribution ^ 
along the full length of the flux tube, giving the total mass m in 
the tube as 


/•SN 

/ p(s)dA(s)ds. 

J ss 


( 2 ) 


Here, dA(s) is the nominal cross-sectional area of the tube, and 
sn and ss are the arc coordinates of the northern (Br > 0) and 
southern (Br < 0) tube footpoints, respectively. Conservation of 
magnetic flux requires that 

|B(s)| dA(s) = |B(sn)| dA(sN) = |B(ss)| dA(ss). (3) 
Likewise, conservation of mass requires that 

m = mt, (4) 


where rh is the rate at which mass is being added to the tube, and 
t is the time elapsed since the magnetosphere last began to refilj^ 
Generalizing equation (31) of jTownsend & Owocki j ( [2005] ), the fill¬ 
ing rate m can be related to the global CAK mass-loss rate M via 


rh = 


M 

AttR. 


/ Br(S]Sl) 

\W^\ 


dA(sN) + 


Br(ss) 

|B(ss)| 


dA{ss)^ 


(5) 


With the arbitrarily chosen timescale t setting the overall mass con¬ 
tained in the magnetosphere, and for given stellar parameters and 
magnetic topology, these equations suffice to calculate the density 
everywhere in the magnetosphere. 


3.3 Applying the ARRM Model to a Ori E 

Applying the ARRM model to a Ori E requires a prescription for 
the magnetic field vector B throughout the magnetosphere. To this 
end, we extrapolate from the surface radial field distribution us¬ 
ing the potential-based approach described by jJardine et al.| ( [200^ . 

^ Recently, [Townsend et al.jjloT^ have presented evidence that the pro- 
cess(es) responsible for limiting the total mass in a magnetosphere may be 
more complex than the episodic centrifugal breakouts envisaged in the ap¬ 
pendix of jTownsend & OwockiH200^ ; however, in the present analysis we 
assume that the latter picture remains correct. 
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Figure 11. Upper plot: The emergent flux from a reference model atmosphere with slightly overabundant chemical composition (shc = —1.0, 6c = —4.1, 
6si = —4.1, and epe = —4.7). Lower plot: The emergent flux from the model atmospheres with modifled abundance of individual elements minus the flux 
from the reference model. A magnifled section of the lower plot showing the redistribution of flux to longer wavelengths is over-plotted in the small upper 
panel. All fluxes were smoothed by a Gaussian Alter with a dispersion of 10 A to show the changes in continuum, a major contribution to SED variability. 



"^ross 


Figure 10. Upper plot: The dependence of temperature on the Rosseland 
optical depth Tross in a reference model atmosphere with a slightly over¬ 
abundant chemical composition (6He = —1.0, 6c = — 4.1, 6si = —4.1, 
and 6Fe = —4.7). Lower plot: The relative change of the temperature in the 
model atmospheres with modifled abundance. 

This method is employed with the assumption that the field remains 
strongly dominant over the wind to a distance within which the ma¬ 
jority of the material is located. As can be seen directly in Fig. 9, 
the material responsible for the Ha emission is located in a radius 


range of ~2-4 R*. Since the density is a strong function of radius, 
we expect the amount of material in the magnetosphere to greatly 
decrease with distance from the star. This agrees with the bulk of 
the material being located close to the Keplerian co-rotation ra¬ 
dius, which is consistent with the results of MHD simulations by 
|ud-Doula et aL| ( [2008] l. Moreover, at this distance, the magnetic ten¬ 
sion dominates the centrifugal force by a factor of ~ 10^, assuring 
minimal distortion from a pure dipole form. The determination of 
where the balance between these forces becomes more equal, and 
how it determines the size and shape of the disc are issues that, 
along with further details about the implementation of this extrap¬ 
olation, will be given in a forthcoming paper (Townsend et ah, in 
prep.). 

Figure [^presents visualizations of the magnetic topology and 
the associated density distribution predicted by the ARRM model 
for a Ori E, at the same five rotation phases illustrated in Figs. 2 
and 4, for an inclination angle i = 75°. The material in the mag¬ 
netosphere is distributed in an oblique, warped disk-like structure, 
qualitatively the same as shown Fig. 4 of [Townsend, Owocki &| 
|Groote| ( [2005| ) for a simple oblique dipole. This is not unexpected 
— the rapid falloff of higher-order components means that the field 
in the ARRM model, at the distances r > 2R^ where the mate¬ 
rial is situated, is close to dipolar. The model simulations keep free 
the parameter which designates the total amount of material in the 
magnetosphere, so that the observations can be best matched. The 
model used here indicates a mass density within the magnetosphere 
of 5 X 10“^^ g cm“^, or a number density of 5 x 10^^ cm“^. 

The light curve and Ha dynamic spectrum predicted by the 
ARRM model are compared against the corresponding observa¬ 
tional data in Figs. and The photometric data are Stromgren 
u band taken by [H^ser ^ al.j ( |1977j ). The observed Ha spectra 
are part of the data presented here and in Paper I. The approach 
used to calculate these observables is similar to that described by 
[Townsend, Owocki & Grootej ( |2005j ), except that we now include 








































10 M.E. Oksala et al. 


Table 1. Ad op ted physical parameters of a Ori E. Table references: ^Oksala| 
|et al.|j20T^ , ^Townsend et al.|j2010} . 


Effective temperature Teff ^ 

Surface gravity log g (cgs)^ 
Projected rotational velocity Urot sin 
Rotational period Prot^ 

Helium abundance 
Carbon abundance 
Silicon abundance 
Iron abundance 


23 000 ±3 000 K 
4.0 ±0.5 
140± lOkms-i 
1^1908229 ± 1^44 x IO-^E; 
-1.1 < £He < 0.6 
-5.0 < £c < -4.0 
—5.9 < 6si < —3.2 
—5.7 < spe < —4.0 


an NLTE photospheric line profile (as described in Sec. 4.1 of Pa¬ 
per I) in the spectral synthesis rather than a flat continuum. The 
comparison has been shown for both i = 75 and 85°; we do not 
show similar plots for z = 55 or 65°, as they make the compari¬ 
son much worse, and as such can be ruled out as possible inclina¬ 
tion angles. Clearly, the model predictions computed with z = 85° 
more closely match the timing of the observed variability, however, 
a clear difference between observations and model can be seen in 
Fig.[^ the secondary minimum at 0 ~ 0.4 is deeper than the pri¬ 
mary minimum at 0 ~ 0.0, the opposite of what is seen in the 
observations. This reversal is also apparent in the strength of the 
Ha emission peaks, with the blue peak being stronger than the red 
one at 0 ~ 0.75 in the observations, but vice-versa for the model. 
We also note that the models do not extend the Ha emission out 
far enough in velocity, particularly on the blue side, something that 
will be explored in the next iteration of this model. Furthermore, 
the improvements added to develop the ARRM model are not suf¬ 
ficient to explain the enhanced emission at phase 0.6. 


4 PHOTOSPHERIC CONTRIBUTION TO THE OPTICAL 
BRIGHTNESS VARIATION 

The previous section (Sec |3.3| ) demonstrates that the current ver¬ 
sion of the RRM model, the ARRM model, is not effective at ex¬ 
plaining the totality of the features seen in the optical variability of 
a Ori E. However, the model is strictly accounting for the circum- 
stellar contribution to these observables, without consideration of 
any other possible contribution to the variability. In the following 
sections, we explore the contribution from the stellar photosphere, 
accounting for the effects of an inhomogeneous abundance distri¬ 
butions of several elements. 


4.1 Simulation of the SED variability 

4.1.1 Model atmospheres and synthetic spectra 


The simulation of the ultraviolet and visual SED variability primar¬ 
ily utilizes the NLTE model atmosphere code TLUS TY (|Hubeny| 
|1988[ [Hubeny & Lanz|p~992| 11995) |Lanz & Hubeny||2003| ). The 


atomic data (adopted from Lanz & Hubeny|2007| l are appropriate 

for B-type stars. These data were originally calculated within the 
Opacity and Iron Projects ( [Seaton et al.|1992[[Hummer et al.|1993| l. 
The calculations assume fixed values of the effective temperature 
and surface gravity for a Ori E (according to Table [^, and adopt a 
generic value of the microturbulent velocity Uturb = 2kms“^. The 
abundances of helium, carbon, silicon, and iron were taken from 
the four-parameter abundance grid given in Table covering the 


Table 2. Individual abundances £He. ^Si. and spe of the model grid 


He 

-1.0 

-0.5 

0.0 

0.5 

1.0 



C 

-5.1 

-4.6 

-4.1 

-3.6 




Si 

-6.1 

-5.6 

-5.1 

-4.6 

-4.1 

-3.6 

-3.1 

Ee 

-5.7 

-5.2 

-4.7 

-4.2 

-3.7 





wavelength X [A] 

Figure 12. The magnitude difference Attza between the emergent fluxes 
calculated with enhanced abundance of individual elements and the refer¬ 
ence flux (see Eq.l^. The fluxes were smoothed by a Gaussian Alter 
with a dispersion of 10()A. 


range of abundances found on the a Ori E surface as inferred from 
the MDI mapping discussed in Sections |23] and [Ad] Abundances 
relative to hydrogen were used, i.e., £ei = log (Aei/An). All other 
elements were set to their solar ( jAsplund et al.|2005] l values. 

Synthetic spectra were calculated using the SYNSPEC code 
for the same parameters (effective temperature, surface gravity, and 
abundances) and transitions as the model atmosphere calculations. 
Additionally, we included lines of all elements with atomic number 
Z ^ 30 not included in the model atmosphere calculation. Angle- 
dependent intensities, /(A, shc, £c, ^si, ^Fe), were computed for 
20 equidistantly spaced values of cos 0, where 0 is the angle be¬ 
tween the normal to the surface and the line of sight. 

The generation of the complete grid of the model atmospheres 
and angle-dependent intensities (Table would require calcula¬ 
tion of 700 model atmospheres and synthetic spectra. Because not 
all abundance combinations correspond to physical conditions de¬ 
termined by the abundance maps, we calculated only those models 
that were realistic. This helped us to reduce the number of calcu¬ 
lated models by about two thirds. 


4.1.2 Calculation of phase dependent SED 

The radiative flux transmitted through a photometric band c at a 
distance D from a star with radius is ( |Mihalas|1978| ) 

/c=(±) J ic{e,n)cosedn. ( 6 ) 

visible 

surface 

The specific band intensity Ic{0,Cl) is obtained by interpolating 
between intensities lc{0, shc, £c, ^si, ^Fe) for abundances found at 
the surface point with spherical coordinates Q. These intensities are 
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Figure 13. Predicted light variations of a Ori E in the Stromgren photomet¬ 
ric system calculated considering only the variations of individual elements. 
The abundance of other elements was fixed. Light curves in individual fil¬ 
ters were vertically shifted to better demonstrate the light variability. 

calculated from the grid of synthetic spectra (see Table|^ according 
to 

POO 

Shq, Sc, Ssi, — / ^c(A)/(A, £He, £c, £Si, ^Fe) dA. 

(7) 

The response function ^c(A) for individual bands is derived by 
either fitting the tabulated response functions or simply assuming a 
Gaussian function. 

The magnitude difference in a given band is defined as 

Arne = -2.5 log , (8) 

where fc is calculated from Eq. and is the reference flux 
obtained under the condition that the mean magnitude difference 
over the rotational period is zero. 

4.2 Influence of abundance on emergent flux 

The bound-free and bound-bound transitions of individual elements 
can modify the temperature of model atmospheres. This can be 
seen in Fig. where we compare the temperature distribution of 
model atmospheres for typical overabundances of a Ori Eto the 
distribution for abundances slightly higher than solar. The bound- 
free (due to ionization of helium, carbon, and silicon) and bound- 
bound (line transition of iron) transitions absorb the stellar radia¬ 
tion, consequently the temperature in the continuum forming re¬ 
gion (xross ~ 0.1 — 1) increases with increasing abundance of these 
elements. In the uppermost layers tross ^ 10“^, the stronger line 
cooling in the models with enhanced abundance leads to the de¬ 
crease of the atmosphere temperature. 

In atmospheres with overabundant helium, carbon, silicon or 
iron, enhanced opacity leads to the redistribution of the flux from 
the short-wavelength part of the UV spectrum to longer wave¬ 
lengths, and also to the visible spectral regions (see Figure [TT] and 



phase 0 

Figure 14. Predicted light variations of a Ori E (solid lines) computed tak¬ 
ing into account helium, silicon, carbon, and iron surface abundance dis¬ 
tributions. Light curves in individual filters were vertically shifted to better 
demonstrate the light variability. The top three light curves are from the 
Stromgren photometric system, with the data again that of jHesser et al.| 
\\911) . The bottom light curve is calculated for the wide-band visual filter 
corresponding to the MOST data jTownsend et al. 12013) . 


[Krticka et al.|2007| l. Consequently, overabundant spots are typically 
bright in the visual and near UV bands, and are dark in the far UV 
bands. 

These flux changes can be detected in the variability of the 
SED and particularly in the variability of visible light. To demon¬ 
strate the optical variations, we plot (Fig.[^ the relative magnitude 
difference 

Am. = -2.5log 

against wavelength. Here is the reference flux calculated for 
slightly overabundant chemical composition (with sne = —1.0, 
£c — —4.1, £si = —4.1, and = —4.7). As can be seen in 
Figure the absolute value of the relative magnitude difference 
is the largest in the u-band of the Stromgren photometric system, 
while it is nearly the same in the visual bands, u, b, and y-band. 
The strong maximum in the u-band model with enhanced helium is 
caused by the filling of the Balmer jump. In the helium rich models, 
the jumps due to hydrogen diminish, while the jumps due to helium 
become visible. 























magnitude [mag] 
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phase (j) 

Figure 15. Comparison of the predicted (solid line) and observed (dots) 
UV light variations for different wavelengths. Curves for individual wave¬ 
lengths were vertically shifted to better demonstrate the variability. 


4.3 Predicted visual light variations 

Predicted light curves are calculated from the surface abundance 
maps derived using MDI in Section 2 and from the emergent fluxes 
computed with the SYNSPEC code, applying Eq.j^for individual 
rotational phases. 

To study the influence of individual elements separately, we 
first calculated the light variations with the abundance map of a 
single element (Eig.[^, assuming a fixed abundance of other ele¬ 
ments (enc = -1.0, £c = -4.1, 5si = -4.1, ^Fe = -4.7). Eig.p^ 
shows that helium, iron, and silicon contribute predominantly to the 
light variations, while the contribution of carbon is only marginal. 
This is due to the large overabundance of these elements within the 
spots and their large abundance variations on the stellar surface. 
The amplitude of the light variations is the largest in the Stromgren 
u-band, as expected from Eig. Because the overabundant re¬ 
gions are brighter in the uvby colors, the predicted light variations 
reflect the equivalent width variations. The light maximum occurs 
at the same phase at which the equivalent width of a given ele¬ 
ment is the greatest. Helium lines are strongest for phase 0.8, con¬ 
sequently the light curve due to helium shows its maximum at phase 
0.8. Iron and silicon lines are strongest for phase 0.6, consequently 
the light curve due to both iron and silicon shows its maximum at 
phase 0.6. 

Incorporating the surface distributions of helium, carbon, sili¬ 
con, and iron in the calculation of light curves (Eig.[^, we obtain 
the results given in Eigure Helium, due to its large overabun¬ 
dance, dominantly influences the light curve of a Ori E. This is the 
opposite behavior seen for HD 37776, where silicon was the main 
cause of the light variability ( [Krticka et al.|2007| ). This is due to the 
higher silicon abundance in HD 37776 by about 0.8 dex. 

The observed visual light curve of a Ori E has a maximum at 
phase ^0.6, when the iron lines are the strongest (Eig[T^, however, 
the model indicates the light curve has a maximum for the phase of 
about 0.8 when the helium lines are the strongest. Consequently, the 
observations strongly disagree with theory in the optical region, and 
the excess brightness at phase 0.6 cannot be explained as an effect 
of inhomogeneous surface abundance of considered elements. 


5 EXTENDED ANALYSIS OF THE UV VARIABILITY 


The light variability caused by the uneven distribution of elements 
on the surface of a Ori E is connected with the flux redistribution 
from the far UV to the near UV and visible regions (see e.g., Eig¬ 
ure [TTJ. Consequently, the light variability in the far UV region pro¬ 
vides an important test of the light variability mechanism. To test 
these predictions quantitatively, we extracted lUE observations of 
a Ori E from the INES databas e ( |Wamstek er^^]J2 000| l using the 
SPLAT package ( Draper 2004 see also Skoda|2008| l. Here ^ 


low-dispersion large aperture spectra in the domain 1250-1900 A 
(SWP camera). 


5.1 Narrow band UV variations 

As a first comparison of the UV fluxes, we focused on narrow band 
variations. The observed and predicted fluxes were both smoothed 
with a Gaussian Alter with a dispersion of 25 A. The resulting pre¬ 
dicted and observed UV variations are given in Pig. Unlike 
the visible region, light curves calculated from the MDI-derived 
abundance maps are able to explain at least the main trends of UV 
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Figure 16. Comparison of observed (lUE) and predicted fluxes in the UV region for selected phases. The predicted flux was smoothed by a Gaussian Alter 
with dispersion 1.3 A, which roughly corresponds the the lUE data. 


light variability of a Ori E. At 1250 A and 1900 A light variabil¬ 
ity is dominated by light absorption in the circumstellar environ¬ 
ment. The observed light curve is nearly constant in the periods 
between the occultations, with the calculated light curve also pro¬ 
ducing little variability. Conversely, there are signiflcant variations 
of the observed flux at 1400 A and 1550 A. These variations can be 
relatively nicely reproduced by our models. The light variability in 
these regions is mainly due to silicon and iron, while helium does 
not signiflcantly influence the light variability in the UV region. 

5.2 UV flux variations 

The variability of the UV spectral energy distribution in Figure 
is only marginal. The predicted and observed emergent flux dis¬ 
tribution agree relatively well during {if = 0.077) and outside 


{if — 0.250, if — 0.836) eclipses. The flux is greatly influenced by 
numerous iron lines, whose accumulation around 1550 A causes 
flux depression in this region. A detailed inspection reveals that 
the predicted flux is somewhat higher than observed in the re¬ 
gion 1150 — 1250 A for all phases. Moreover, the predicted flux 
is slightly lower than observed in the region 1650 — 1900 A during 
the maximum of the UV flux variability {f — 0.836). This may 
possibly be the source of disagreement between the light curves in 
the optical region. 

5.3 UV line variations 

The shapes of many weak lines and the resonance lines of Si IV at 
~ 1400 A seem to be well reproduced by the models in Figure [T^ 
On the other hand, the strengths of strong lines of Si II and Si III 
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at 1270 A and 1300 A and the C IV resonance lines are not repro- 
duced by the models. We use the following indices to study line 


variability: 


a(Cll 1336 A) = ?ni335.2 — - (mi332.o + ^isss.rs) , 

(lOa) 

a(ClII 1175 A) = ?nii75.7 — - (mneg.o + ^ 1181 . 25 ) , 

(10b) 

a(ClII 1247 A) = 7711247 .5 — 2 (^1244.45 + 7771249.15) 

, (10c) 

a(ClV 1548 A) = 7771547.9 — - (7771529.45 + 7771557.95) 

, (lOd) 

a(SilV) = 7771393.65 “ ^ (^1388.65 + 7771398. 7 ) 

. (lOe) 

These indices were used in the analysis of|Krticka et al.|(|2013|) 
for HD 64740, who modified a similar set of indices developed by 
[Shore, Brown & Sonnebornj(|1987j). Eigure pTj compares the ob- 


served and predicted variations of the line indices with respect to 
their mean values. We are able to explain the variations of many 
individual carbon and silicon lines. The disagreement between pre¬ 
dicted and observed line profiles of Si ll and Si III could be caused 
by an inhomogeneous vertical distribution of silicon in the atmo¬ 
sphere. The observed line strengths of C IV at 1548 A and 1551 A 
are significantly larger than predicted. The observed asymmetry of 
these lines with extended blue wings at each phase indicates forma¬ 
tion in the wind. The line strength variations result from the vari¬ 
ations of wind mass-flux on the stellar surface. Enhanced absorp¬ 
tion in the C II line at 1336 A, C IV lines, and Si IV lines observed 
around the phase (/? 0 possibly originates in the corotating mag- 

netospheric clouds. 


6 CONFRONTATION: ROUND 2 

In paper I, we first confronted the RRM model with new high- 
resolution spectropolarimetry of a Ori E to demonstrate that 
the field configuration assumed by [Townsend, Owocki & Groot^ 
( |2005| ) was incompatible with current observations. This first step 
revealed that more comprehensive modeling was required to ex¬ 
plain the intricate physical phenomena of this star. While the initial 
success of that model is now attributed to a bug in the code, its in¬ 
spiration has driven studies to better understand the intimate details 
of the interaction of rotation, magnetic field, and mass-loss, pro¬ 
pelling the new wave of high quality data acquisition presented in 
this and other papers. The current version of the RRM model now 
does a marginally good job of replicating the qualitative variations 
of the system, but several discrepancies appear to test the accuracy 
of the model, and the excess brightness in the optical photometry 
at phase 0.6 persists. 

The magnetic field topology in this application of the ARRM 
model is taken from the MDI analysis in Section [2^ and is there¬ 
fore consistent with the observations of longitudinal magnetic field. 
Several checks were performed to to ensure that the model adopted 
the same phasing and rotation as the derived maps, particularly af¬ 
ter the discovery of the bug in the initial model. The Ha emission 
comparison is shown in Eigure [9] and while the model reproduces 
much of the variability, there is one glaring inconsistency in the 
strength of the red-shifted cloud at phase ~0.75. If we carefully in¬ 
spect the other features throughout the rotation, it seems this one as¬ 
pect is the only serious problem, indicating the model is predicting 
too much emitting material in that “cloud” at that particular phase. 
Likewise, Eigure shows a similar reversal in strength of the the 
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Figure 17. Observed variability of the line indices (dots) in comparison 
with the predicted indices (lines). Upper panel: Carbon line indices. Lower 
plot: Silicon line index. The line indices were calculated according to Eq. 
(10). 


model photometric minima with the deeper minimum occurring at 
the opposite phase as compared to the observed light curve. 

Section [43] demonstrates clearly there is no photometric reso¬ 
lution to the issue of the excess brightening in the second half of the 
rotation phase, or the irregular shape of the minimum at phase 0.4, 
both clearly seen in the MOST photometric light curve. As model¬ 
ing the photometric contribution to the brightness variability does 
not resolve the discrepancies, we must at this point consider that 
the ARRM model is not currently capable of explaining all of the 
observed variabilities of a Ori E and similar stars in detail, and that 
additional physics is necessary. 


6.1 Where is the missing physics? 

Where does this “extra” brightness come from? If we look at Eig¬ 
ure geometry could be a possible explanation. Near phase 0.6- 
0.8, where we see brightening, the magnetosphere is seen more 
or less pole-on, at which we are able to view all of the emitting 
surface of the magnetospheric material. At the phase 0.2, we see 
that the material is viewed more edge-on, with a much less direct 
view of most of the emitting surface. It is possible (and likely) that 
the star’s radiation may be reflecting off of the inner edge of the 
magnetosphere, thus creating a brightening effect. At this iteration, 
this and other types of scattering are not considered in the RRM 
model, but likely contribute to the observed brightness. This is ex¬ 
tra physics that would need to be implemented into the next version 
of the code to diagnose the importance and type of scattering that 
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could be responsible, i.e., Rayleigh scattering vs. electron scatter¬ 
ing. 

Another consideration is the shape of the “clouds”. The 
ARRM model predicts, for a Ori E, two higher density concen¬ 
trations within a larger disk-like structure. However, the broadband 
polarization study of |Carciofi et al.|f2013| ) suggests that the ma¬ 
terial is configured in a much more concentrated structure, much 
more similar to dumbbells that a disk. A more physically realistic 
density distribution may present different model predictions. An¬ 
other effect of the interaction of the star and the magnetic field is the 
way in which material leaves the surface of the star. In a “normal” 
OB star, the mass-loss is purely dependent on the strength of the 
wind, and is proportional to the luminosity. Although the ARRM 
model considers effects from the magnetic field, the model still as¬ 
sumes standard CAK relations. The influence of strong field lines 
may modify the mass-loss process, changing the rate at which the 
material leaves the star resulting in higher or lower mass-loss. The 
complexity of the field structure may further constrict movement of 
mass from the surface. An effect like this could be responsible for 
the discrepancies between the observed Ha emission and the model 
prediction. Similarly, the influence of the surface abundance inho¬ 
mogeneities may play an additional role in altering the mass-loss 
rate. 


7 SUMMARY AND FUTURE WORK 

In the first paper of this series, we evaluated the appropriateness of 
the magnetic field assumed by fTownsend, Owocki & Groote| ( |2005| > 
in which they applied the RRM model to the physical case of a Ori 
E. The conclusion of that work, that a more complex field was re¬ 
quired, served as the motivation for the analyses contained in this 
subsequent paper. We have used the same high-resolution spectra 
from Paper I to produce MDI maps of both the magnetic topology 
and the surface abundance distributions of He, C, Si, and Ee. The 
derived magnetic topology, which can be represented by a super¬ 
position of a dipole plus quadrupole, is consistent with our previ¬ 
ous assertion that the field contains higher-order components. This 
magnetic configuration was then fed into the newly developed ar¬ 
bitrary magnetic field version of the RRM model to re-evaluate the 
consistency between the model predictions and the observed varia¬ 
tions. The resultant predictions qualitatively resemble the observed 
spectral and photometric variations, but fail to reproduce still sev¬ 
eral details. While the magnetic field is fixed, the model now pre¬ 
dicts a larger density in one of the clouds, which is observed, but 
located on the wrong side of the star. The new model also does 
not resolve the previously-identified issues, specifically the excess 
brightness observed at phase ~0.6. 

While clearly some of these issues are related to failings of the 
model, we proceeded to use the second result of the MDI modeling, 
the abundance maps, to determine the effect of the surface inhomo¬ 
geneities on the photometric light curve. With these maps as input, 
a synthetic light curve was computed with the aid of a grid of stellar 
atmosphere models. The total resultant light curve, a product of the 
contributions from He, C, Si, and Ee, provides no solution for the 
discrepancies between the model and the observations. The light 
curve is dominated by the He contribution, as this element exhibits 
the largest variation of abundance across the surface, however, its 
maximum abundance is viewed at phase 0.8, too late to be helpful 
in accounting for the feature at phase 0.6. An analysis of the UV 
light variability shows good agreement between observed variabil¬ 
ity and computed light curves, however, some narrow wavelength 


bands are affected by the eclipsing of the magnetospheric material, 
causing enhanced absorption at similar phases as in the optical. 

We thus conclude that the discrepancy between the models and 
the observations lies in the treatment of the stellar magnetosphere, 
and that the RRM model must undergo further revision. Adding 
the ability to accommodate an arbitrary magnetic field topology is 
just the first step towards our greater understanding of this system, 
but there are multiple issues to be addressed in future iterations. 
The inclination angle first and foremost should be determined using 
all the available information. The model will need to incorporate 
and test the effects of scattering processes so as to determine their 
ability to alter the synthetic observations. Tests should be run to 
understand how the field interacting with the wind could change the 
stellar mass-loss. When the RRM model was developed, massive 
star magnetism was just beginning to progress, and at the time, a 
Ori E was a unique object, the only known member of its class. 
However, in the past decade, research has grown rapidly, and there 
are a growing number of stars discovered to have similar properties, 
and exhibit similar variabilities. By optimizing the RRM model for 
cr Ori E, we can create a blueprint to fully comprehend the physical 
properties, the interactions within and outside of the star, and the 
evolutionary implications for this fascinating group of stars. 
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